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Abstract: The model of a tired random walker, whose jump-length decays 
exponentially in time, is proposed and the motion of such a tired random 
walker is studied systematically in one, two and three dimensional contin¬ 
uum. In all cases, the diffusive nature of walker, breaks down due to tiring 
which is quite obvious. In one dimension, the distribution of the displace¬ 
ment of a tired walker remains Gaussian (as observed in normal walker) with 
reduced width. In two and three dimensions, the probability distribution 
of displacement becomes nonmonotonic and unimodal. The most probable 
displacement and the deviation reduces as the tiring factor increases. The 
probability of return of a tired walker, decreases as the tiring factor increases 
in one and two dimensions. However, in three dimensions, it is found that 
the probability of return almost insensitive to the tiring factor. The proba¬ 
bility distributions of first return time of a tired random walker do not show 
the scale invariance as observed for a normal walker in continuum. The 
exponents, of such power law distributions of first return time, in all three 
dimensions are estimated for normal walker. The exit probability and the 
probability distribution of hrst passage time are found in all three dimen¬ 
sions. A few results are compared with available analytical calculations for 
normal walker. 
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I. Introduction: 

In statistical physics, process of polymerization[Tl [2], diffusion [3] in re¬ 
stricted geometry etc are some classic phenomena, which have drawn much 
attention of the researcher in last few decades. The underlying mechanism 
of such physical phenomena are tried to explain by random walk [3]. Differ¬ 
ent types of random walk are studied on the lattice in different dimensions 
by the method of computer simulation. The absorbing phase transition in a 
conserved lattice gas with random neighbour particle hopping is studied[5j. 
Quenched averages for self avoiding walks on random lattices [U], asymptotic 
shape of the region visited by an Eulerian walker [7], linear and branched 
avalanches are studied in self avoiding random walks [8], effects of quenching 
is studied in quantum random walk recently [9]. The drift and trapping in 
biased diffusion on disordered lattices is also studied [TO]. 

Very recntly, some more interesting results on random walk were re¬ 
ported. The average number of distinct sites visited by a random walker on 
the random graph[TT], statistics of first passage time of the Browian motion, 
conditioned by maximum value of area|12j are studied recently. It may be 
mentioned here that the first passage time in complex scale invariant media 
was studied [13]. The thoery of mean hrst passage time for jump processes 
are developed m and verified by applying in Levy flights and fractional 
Brownian motion. The statistics of the gap and time interval between the 
highest positions of a Markovian one dimensional random walker [15], the 
universal statistics of longest lasting records of random walks and Levy 
flights are also studied [T6]. 

The random walks in continuum are studied to model real life prob¬ 
lems. The exact solution of a Brownian inchworm model and self-propulsion 
was also studied [T7], theory of continuum random walks and application in 
chemotaxis was developed [18]. Random walks in continuum was also stud¬ 
ied for diffusion and reaction in catalyst [19] . Very recently, the random walk 
in continuum is studied with uniformly distributed random size of flight [20]. 
The statistics of Pearson walk is studied [711122] in two dimensions for shrink¬ 
ing stepsize and found a transition of the endpoint distribution by varying 
the initial stepsize. 

The living random walker in continuum gradually becomes tired as the 
time passes, in reality. This would reduce its energy, as a result the size of 
flight gets reduced gradually with time. The first-passage properties [21] of 
a walker is important in various aspects, namely, the fluorescence quenching 
in which a fluoresecent molecule stops while reacting with a quencher, firing 
neurons when the fluctuating voltage level first reaches a specified value, 
in econophysics, the execution of buy/sale orders when a stock price first 
reaches a threshold. What will be the first passage properties if the stepsize 
of a Pearson walker decreases exponentially in time ? In this paper, address¬ 
ing this particular problem, a model of tired random walker is proposed in 
continuum and statistics of its motion is studied systematically in one, two 
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and three dimensional continuum. The first passage properties, return and 
exit probabilities are studied here. The numerical results of detailed statis¬ 
tical analysis of the motion of a tired random walker is also reported here. 
This paper is organised as follows: In the next section (section-II) the model 
of tired random walk is proposed and the results obtained from numerical 
simulations are given. The paper ends with a summary given in section-III. 
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II. Model and Results: 

Generally, the motion of a random walker is studied by considering the 
time (t) independent size (i?) of flight in each move. In this study, the model 
of a tired random walker is proposed in such a way that the size of flight of 
a walker decreases exponentially as R{t) = e“"*. A simple logic behind it 
may be stated as follows; if a living cell is moving continuously, its energy 
(basically kinetic energy) gradually decreases and hence the velocity, which 
in turn reduces its size of flight (i.e., jump-length per unit time). Here, a 
is tiring factor. The statistical behaviour of such a tired random walker is 
studied in one, two and three dimensional continuum. It may be noted here 
that such kind of behaviour of a tired random walker cannot be studied on 
the lattice. 

In one dimension, the size of flight in each time step is xi{t) = e“"*. 
A walker starts its journey from the origin having the equal probability of 
choosing the left and right direction. The updating rule, in one dimemsional 
tired walk, may be expressed as: x{t -|- 1) = x{t) ± xi{t) 

In two dimensions (planar continuum), the tired walker starts its journey 
from the origin and it has a uniform probability of choosing any random 
direction (6) distributed between 0 and 2tt. Its motion can be represented 
mathematically as: 

x{t -|- 1) = x{t) + R{t)cos{6)] y{t -|- 1) = y{t) + R{t)sin{9) (1) 

The displacement at time t is r{t) = x^{t) + y‘^{t). In planar continuum, 

the area of the region visited by a tired walker is obviously shorter than 
that visted by a normal walker, in a specified course of time. A typical 
such comparison is shown in Fig-1 with a = 0.001. As a result, the mean 
square displacement does not show diffusive behaviour as shown by a normal 
walker. In long time, it gets saturated (motion stops practically). 

A typical such comparison is shown in Fig-2 for a = 0.001 and a = 
0.0005. The similar behaviours are also observed in one and three dimensions 
(not shown). The tired walk is not diffusive (< >~ t) as observed in 

normal (a = 0) walk. It is also observed that the motion stops earlier if the 
tiring factor a increases. 

Now, let it be discussed systematically in one, two and three dimen¬ 
sions. In one dimensions, the probability distribution of the displacements 
of a walker are studied for a = 0.0 (normal), 0.001 (moderately tired) and 
0.01(heavily tired). As usual, the distribution is normal (Gaussian) with 
zero mean in all the cases. However, as the tiring factor increases the distri¬ 
bution becomes sharper and sharper. These are depicted in Fig-3. Here, it 
may be mentioned that the values of a and the maximum time allowed (Nt) 
are such that the walker gets frozen (due to exponential decrease of step- 
size after such long time). The distribution shown in Fig-3, is practically 
the density distribution of frozen walker. It would be interesting to study 
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the density distribution of these frozen walker as a function of a through 
the scaling. 

What will be the probability of return(PR) in one dimension ? First of 
all, in continuum one should be careful in defining the probability of return. 
In the lattice the probability of return is defined as the walker returns to its 
initial starting point. However, in continuum, it is quite unlikely that a tired 
walker returns to its initial starting point. Here, one may think that whether 
the tired walker returns within a linear zone {[—rz,rz\) centered around the 
origin. Now put a large number (W) of walker at origin and allow them to 
walk (with different random sequence) upto a certain time (W) and then 
check how many walkers return within the preassigned returning zone (of 
size rz)- The calculated fraction is the probability of return (within time of 
observation Nt) in this particular model. In the lattice model this probability 
is 1, which can also be derived from exact calculations [23| . In this model of 
tired walker, considering = 0.5, Pr = 0.992 for a = 0.0. This numerical 
estimate of returm probability agrees well with exact calculation of return 
probability {Pr = 1)[23] in one dimensional normal random walk. It may 
be noted here that for a = 0, the walker returns at origin (the starting 
point also) and the probability of return can be compared to that obtained 
in random walk on one dimensional lattice. As the tiring factor increases, 
Pr decreases. For moderately tired (a = 0.001) walker, Pr = 0.955 and for 
heavily tired walker (a = 0.01), Pr = 0.874. In Fig-4, the Pr is plotted 
against Nt for various values of a. Now, this probability of return (Pr) must 
depend on the size (rz) of returning zone. To study the dependences of Pr 
on rz, Pr is studied as a function of for different values of a and shown 
in Fig-5. It shows that the Pr grows as rz increases in the case of tired 
walker (a = 0.001,0.01), but Pr does not depend on for normal walker. 
It is important to note here that even for heavily tired (a = 0.01) random 
walker, the size of the flight, after t = 10 is larger than 0.90. So, the range of 
values of rz, chosen here, does not have any chance that the walker remains 
in the returning zone immediately after starting its journey. So, the choice 
rz =0.5 is quite safe to study the probability of return in this context. 

How long a tired walker takes to return first time within returing zone 
? How does the distribution of this first returning time (t^) look like ? The 
probability distribution of first returning time (tr) of a tired random walker is 
shown in Fig-6. A normal walker (a = 0.0) shows a scale invariant {P{tr) ~ 
tr^) distribution of first returning time (tr)- The exponent estimated is 
(3 ~ 1.49. This result agrees well with analytical result [Mj, where it is 
found P(tr) ~ tr ' • However, this scale invariant nature of the distribution 
of first returning time, breaks down in the cases of tired walking (for a = 
0.001,0.01) (see Fig-6). More detail investigation is required to propose any 
functional behaviour of P{tr) for a ^ 0. 

In one dimension, how long (tg) a tired walker takes to exit (first time) 
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from a zone ([—rg, re]) ? The probability distribution [P{te)) of first passage 
time (tg) (for a fixed value of rg = 25.0), is studied for different values of a 
and shown in Fig-7(a). As the a increases, the most probable first passage 
time decreases. It should be noted here that the probability of first passage 
is defined (in this study) as the probability to escape in a given time, from a 
bounded ([—rg,rg]) linear (in one dimension) region. If it would be defined 
as the probability to escape through a given point (say x = re) the power 
law {P{te) ~ distribution in long time limit (tg —>■ oo) is found which 

supports the analytical prediction {P{te) ~ tg^^^)[2l]- This is shown in 
Fig-7(b). 

What is the probability of exit (Pg) of a tired walker one dimensional 
continuum ? The exit probability (for a fixed time of observation Nt) from a 
zone of absolute distance rg (measured from the origin) is also studied, in one 
dimension, as a function of rg and shown in Fig-8. Here, the exit probability, 
of a tired walker, was found to decreases as the absolute distance of zone 
(rg) increases. However, it remains fixed (nearly 1) for a normal walker [25]. 
It may also be noted that, the rate of fall of exit probability increases as the 
tiring factor (a) increases. 

In two dimensions, the motion of a tired random walker is studied by us¬ 
ing the rule given in equation -1. Here, the mean square displacement < r^ > 
is proportinal to the time t for a = 0.0, reveals the conventional diffusive 
(< r^ >~ t) behaviour. However, a moderately tired (a = O.OOl&O.OOOS) 
walker does not show long time diffusive behaviour. This is quite obvious 
and already shown in Fig-2. 

The distribution of absolute displacement is nonmonotonic unimodal 
function. It is shown in Fig-9. It is observed that the maximum probability 
of finding the walker at a distance (r^) from the origin and the average 
distance (f) both decreases as the tiring factor (a) increases. In this case, 
Xm = 64.0, f = 89.13 for a = 0.0, = 15.0, f = 20.04 for a = 0.001 and 

rm = 5.0, f = 6.42 for a = 0.01. Here also this distribution is practically 
the density distribution of frozen walker. 

What will be the probability of return {Pr) of a tired walker in planar 
continuum ? The probability of return within a circle of return having radius 
rz = 0.5 is studied as a function of maximum time of observation Nt and 
shown in Fig-10. In planar continuum, a tired walker has a probability of 
return in a circle of radius = 0.5 as follows: for a = 0.0, Pr = 0.737, for 
a = 0.001, Pr = 0.620 and for a = 0.01, Pr = 0.524. 

For a fixed value of Nt, the probability of return of a tired walker in 
planar contunuum, grows as the radius of returning zone increases. This is 
shown in Fig-11. 

Here, like in one dimensional tired walker, the probability distribu¬ 
tion of first returning time (p) shows a scale invariance {P{tr) ~ tr^) 
for a = 0.0 with [3 ~ 1.10. However, the analytic result[23] suggests 
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P{tr) ~ l/{tr{ln{tr)‘^). The possible reason of disagreement may be stated 
as follows: in the analytic calculation of P{tr), it was defined as the proba¬ 
bility of return exactly at the origin from where the walker has started its 
journey. However, in the numerical simulation, P{tr) is defined as the prob¬ 
ability of return (first time) within a circular zone of radius r^. As the tiring 
factor (a 7 ^= 0 ) increases, the scale invariance nature of the distribution on 
first returning time breaks down. This is demonstrated in Fig-12. 

In two dimensions, the distribution of first passage time (for a fixed 
distance rg = 25.0), is studied for different values of a and shown in Fig- 
13. As the a increases, the most probable first passage time and mean first 
passage time decreases. 

The exit probability (for a fixed time of observation Nt) from a circular 
zone of radius Vg (measured from the origin) is also studied, in two dimen¬ 
sions, as a function of Vg and shown in Fig-14. Here, the exit probability was 
found to decreases as the radius of circular zone (vg) increases. Here also, 
the rate of fall of exit probability increases as the tiring factor (a) increases. 
However, the exit probability of a normal walker [a = 0) remains unchanged 
(nearly 1 ) as rg increases. 

The tired walk in three dimensional continuum can also be generalized. 
The updating of coordinates obey the following rule; 

x{t -|- 1) = x{t) + R{t)sin{9)cos{(j))-, y{t -|- 1) = y{t) + R{t)sin{6)sin{4>) 

z{t -|- 1) = z{t) + R{t)cos{6) (2) 

where R{t) = 9 is uniformly distributed random angle between 0 and 

TT and (j) is uniformly distributed random angle between 0 and 27r. The 
displacement at time t is r{t) = ^Jx‘^{t) + y‘^{t) + z‘^{t). 

In 3D continuum, the motion of a tired random walker is studied by 
using the rule given in equation -2. Here, the mean square displacement 
< > is proportinal to time t for a = 0 . 0 , reveals the diffusive behaviour 

(not shown). However, a tired (a 7 ^ 0) walker does not show long time 
diffusive behaviour (not shown). 

The probability distribution of absolute displacement (or the density 
distribution of frozen walker in reality) in 3D continuum is observed to be 
a nonmonotonic unimodal function. It is shown in Fig-15. It is observed 
that the maximum probability of finding the walker at a distance (r^) from 
the origin and the mean displacement (f) both decreases as the tiring factor 
(a) increases. In this case, = 72.0, f = 91.59 for a = 0.0, = 16.0, 

f = 20.65 for a = 0.001 and Vm = 3.0, f = 6.96 for a = 0.01. 

What will be the probability of return in 3D continuum ? The probability 
of return within a sphere of return having radius = 0.5 is studied as a 
function of maximum time of observation Nt and shown in Fig-16. In 3D 
continuum, unlike the cases in ID and 2D continuum, a tired walker has 
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a probability of return in a sphere of radius = 0.5 is almost insensitive 
{Pr = 0.226) of the tiring factor a. 

For a fixed value of Nt, the probability of return of a tired walker in 
3D contunuum, grows as the radius of returning zone increases keeping the 
independence on tiring factor a. This is shown in Fig-17. 

In 3D continuum, the probability distribution of first returning time {tr) 
shows a scale invariance {P{tr) ~ tr^) for a = 0.0. The exponent estmated 
/? ~ 1.48. Accidentally, this is close the analytical prediction {P{tr) ~ 
tr ^'^^)[23]. As the tiring factor {a ^ 0) increases, the scale invariance of the 
distribution on first returning time, breaks down. This is demonstrated in 
Fig-18. 

In three dimensions, the distribution of first passage time (for a fixed 
distance rg = 25.0), is studied for different values of a and shown in Fig-19. 
As the a increases, the most probable first passage time and the mean first 
passage time decreases. 

The exit probability (for a fixed time of observation A*) from a spherical 
zone of radius rg (measured from the origin) is also studied, in two dimen¬ 
sions, as a function of rg and shown in Fig-20. Here, the exit probability, 
of a tired walker, was found to decreases as the radius of circular zone (rg) 
increases. However, like the earlier cases, it reamins same (nearly 1) for 
all Tg. Here also the rate, of fall of the exit probability of a tired walker, 
increases as the tiring factor a increases. 


III. Summary: 

In this article, a model of tired random walker in continuum is proposed. 
Generally, a random walker moves with constant size of flight. However, as 
the time passes, if the walker gets tired, one should think of a time dependent 
size of flight. Here, this size of flight decays exponentially with time. The 
motion of such a tired walker is studied in one, two and three dimensional 
continuum. In this statistical investigation, the distribution of the absolute 
displacement, mean displacement, probability of return (within a specified 
zone), distribution of time of first return are studied systematically. In 
one and two dimensional continuum, the probability of return decreases as 
the tiring factor increases. However, in three dimensional continuum, this 
probability of return seems to be independent of the tiring factor. The 
distribution of first returning time in all dimensions (for normal walker with 
tiring factor a = 0), shows power law behaviours. This scale invariance 
of the distribution of first returning time breaks down for a 7 ^ 0 in all 
dimensions. In the study of first returning probability, a very important 
point should be mentioned. For a = 0, the probability of return could be 
compared with that calculated analytically |24j in one dimension only, where 
the walker can return to the initial point. In higher dimensions, it returns 
within a circular (spherical) zone in two (three) dimensions. 

The exit probability and the distribution of first passage time are studied. 
In all dimensions, the exit probability was found to decrease as the size of 
the zone (from where the tired walker exits out) increases. The rate of 
decrease of the exit probability was found to increase as the tiring factor a 
increases. Here also the probability of first passage (for a = 0) can only be 
compared with analytical calculations |24j in one dimension, if it is defined 
as the probability of escape through a particular point. 

The first passage time is defined (in this simulational study) as the time 
required by a walker to exit from a specified zone. This time has a distribu¬ 
tion and this distribution is studied for various values of a. It was observed, 
in all dimensions, the most probable first passage time decreases as a in¬ 
creases. A rigorous analysis and possible scaling behaviour (if any) may be 
investigated. 

Some more interesting studies can be done in this field. In this paper, 
only the numerical results are reported. A rigorous mathematical formula¬ 
tions of first passage properties for tired walk has to be developed following 
the same already developed |24j for normal walk (a = 0 ). 

The possibilities of scaling of distribution of return time, distribution 
of first passage time, distribution of distances and exit probabilities with 
respect to the tiring factor (a) has also to be explored. 
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Fig-1. A sample random walk in 2D continuum. The top one shows a walk 
for a = 0.0 and the bottom one shows a tired walk for a = 0.001. Here, in 
both figures same sequence of random numbers are used. Note the scales of 
the axes of two figures. Nf = 10^ in both cases. 
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Fig-2. The mean square displacement (< >) versus time (t) in two 

dimensions for various values of a (marked by different symbols). The tired 
walkers {a ^ 0) do not show the normal (< >~ t) diffusive behaviour. 
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Fig-3. The distribution (P(x)) of the displacements (x) of a tired random 
walker in one dimension, (o) represents a = 0.0, (□) represents a = 0.001, 
and (*) represents a = 0.01. Note that tiring reduces the width of the 
distribution. 
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Fig-4. Probability of return {Pr) plotted against the maximum time Nt 
in one dimension. Different symbols correspond to different values of a. 
a = O.O(o), a = 0.001(«) and a = 0.01(*). Here in all cases Ng = 10^. The 
absolute distance (r^) of returning zone ([—is = 0.5 here. 
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Fig-5. Probability of return (Pr) plotted against the absolute distance 
of returning zone {[—Vz, r^]) in one dimension. Different symbols correspond 
to different values of a. a = O.O(o), a = 0.001(«) and a = 0.01(*). 
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Fig-6. Probability distribution {P{tr)) of first returning time (t^) in one 
dimension. Different symbols correspond to different values of a. a = O.O(o), 
a = 0.001(«) and a = 0.01(*). The solid line is y = 
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Fig-7. Probability distribution (P(te))of first passage time (te) in one di¬ 
mension. (a) Different symbols correspond to different values of a. a = 
O.O(o), a = 0.001 (•) and a = 0.004(*). Here, Ng = 10®, = 10^ and 

Tf. = 25.0. In (a) the first passage is defined as the probability of exit first 
from a bounded i[—re,re]) linear region around the origin, (b) the first pas¬ 
sage is defined to cross first a point (here rg = -1-25.0). The solid line is 
y = 1000x“^'‘^® supporting analytical prediction {P{te) ~ 
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Fig-8. Exit probability Pg plotted against Vg in one dimension. Different 
symbols correspond to different values of a. a = O.O(o), a = 0.001(») and 
a = 0.01(*). Here, Ng = 10^ and Nt = 10^. 
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Fig-9. The distribution of displacements (r) of a tired random walker in 
two dimensions, (o) represents a = 0.0, = 64.0, f=89.13, (□) represents 

a = 0.001, Tm = 15.0, f=20.04 and (*) represents a = 0.01, = 5.0, 

f=6.42. 
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Fig-10. Probability of return {Pr) plotted against Nt in two dimensions. 
Different symbols correspond to different values of a. a = O.O(o), a = 
O.OOl(n) and a = 0.01(=i<). Here in all cases Ng = 10®. The radius of 
circular returning zone is = 0.5. 
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Fig-11. Probability of return (Fr) plotted against the radius of returning 
zone {tz) in two dimensions. Different symbols correspond to different values 
of a. a = O.O(o), a = O.OOl(n) and a = 0.01(*). 
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Fig-12. Probability distribution {P{tr))oi first returning time (t^) of two 
dimensions. Different symbols correspond to different values of a. a = 
O.O(o), a = O.OOl(n) and a = 0.01(*). Solid line represents y = 

The dotted line is y = 200/(a;(/n(x))^)[2l]. 
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Fig-13. Probability distribution {P{te)) of first passage time (te) in two 
dimensions. Different symbols correspond to different values of a. a = 
O.O(o), a = 0.001(«) and a = 0.002(*). Here, Ng = 10®, Nt = 10^ and 
re = 25.0. 
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Fig-14. Exit probability Pg plotted against Vg in two dimensions. Different 
symbols correspond to different values of a. a = O.O(o), a = 0.001(») and 
a = 0.01(*). Here, Ng = 10^ and Nt = 10^. 
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Fig-15. The distribution of the displacements of 
in three dimensions, (o) represents a = 0.0, 
represents a = 0.001, = 16.0, f=20.65 and 

Tm = 3.0, f=6.96. 
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Fig-16. Probability of return {Pr) versus Nt in three dimensions. Different 
symbols correspond to different values of a. a = O.O(o), a = O.OOl(n) and 
a = 0.01(*). Here in all cases Ng = 10®. The radius of spherical returning 
zone is = 0.5. 
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Fig-17. Probability of return (Fr) versus in three dimensions. Different 
symbols correspond to different values of a. a = O.O(o), a = O.OOl(n) and 
a = 0 . 01 (*). 
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Fig-18. Probability distribution {P{tr)) of first returning time (tr) in three 
dimensions. Different symbols correspond to different values of a. a = 
O.O(o), a = O.OOl(n) and a = 0.01(*). Solid line represents y = 
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Fig-19. Probability distribution {P{te)) of first passage time (tg) in 
dimensions. Different symbols correspond to different values of a. 
O.O(o), a = 0.001(«) and a = 0.002(*). Here, Ns = 10®, Nt = 10' 
rg = 25.0. 
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Fig-20. Exit probability Pg plotted against Vg in three dimensions. Different 
symbols correspond to different values of a. a = O.O(o), a = 0.001(») and 
a = 0.01(*). Here, Ng = 10^ and Nt = 10^. 
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